Triplet supercurrent in ferromagnetic Josephson junctions by spin injection 
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We show that injecting nonequilibrium spins into the superconducting leads strongly enhances 
the stationary Josephson current through a superconductor-ferromagnet-superconductor junction. 
The resulting long-range super-current through a ferromagnet is carried by triplet Cooper pairs 
that are formed in s-wave superconductors by the combined effects of spin injection and exchange 
interaction. We quantify the exchange interaction in terms of Landau Fermi-liquid factors. The 
magnitude and direction of the long-range Josephson current can be manipulated by varying the 
angles of the injected polarizations with respect to the magnetization in the ferromagnet. 
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INTRODUCTION 



Studies of hybrid structures combining supercon- 
ducting and ferromagnetic components attract much 
attention due to their unique, rich, and complex 
physical properties that are promising in a number 
of potential applications^. The interface of an s- 
wave superconductor with a ferromagnet is charac- 
terized by an unusual proximity effect that is spa- 
tially oscillating and can lead to a sign reversal of 
the critical current through supcrconductor-fcrromagnet- 
superconductor (SFS) Josephson junctions. Such a rever- 
sal is equivalent to a 7r-shift in the current-phase relation 
for the Josephson current. This interesting property is 
a motivation for using the so-called 7r-junctions as el- 
ements of superconducting quantum circuits for poten- 
tial application in quantum computing 2 -. However, the 
proximity effect in ferromagnets docs not reach far. Two 
critical tasks are to extend its range and to find a way to 
manipulate the 7r-junction in order to switch the device 
between its various phase states. In contrast, Cooper- 
pairs can be transferred over relatively long distances 
even in ferromagnets, if they are in a triplet state with 
±1 projections of their total spin onto the spin quantiza- 
tion axis. Various mechanisms have been proposed that 
convert a singlet pair into a triplet pair, such as a spa- 
tially dependent magnetization^, spin-flip scattering at 
FS interfaces^, and precessing magnetization^. A num- 
ber of works in this direction has been reviewed in RefH 

In this work, we will show that these tasks can be ful- 
filled via the production and manipulation of a long-range 
proximity effect by injecting spins into superconducting 
leads. The novelty of our idea is based on the impor- 
tant, and so far unaddrcsscd, role played by the electron- 
electron interaction in SFS. Our insight is that the com- 
bined effects of spin-injection and electron-electron inter- 
action generate a long-range proximity effect despite the 
strong exchange field in the ferromagnet. The conven- 
tional wisdom is that spin polarized electrons can only 
exist as excitations in s-wave superconductors, since the 
Cooper pairs do not carry a spin. However, we will 
demonstrate that this simple picture, which is based on 



the neglect of electron-electron interactions beyond su- 
perconducting pairing correlations, misses qualitatively 
important effects. Quantitatively, in simple metals, the 
exchange interaction of itinerant carriers is noticeable 
and can be described in terms of Landau Fermi- liquid fac- 
tors. Although the exchange interaction docs not cause 
ferromagnetism in s-wave superconductors, it causes a 
transfer of spin polarization from the quasi-particle exci- 
tations to the condensate, in the form of polarized triplet 
Cooper pairs. When such a triplet pairing is generated 
by the combined effects of spin-injection and exchange 
interaction, these pairs subsequently tunnel through the 
ferromagnetic layer via the long-range proximity effect, if 
the spin polarizations in the leads and the layer are not 
collinear. Only at this stage, which includes the so far 
unaddressed important electron-electron interaction, the 
situation becomes similar to proposals of Ref|l||3| where 
an inhomogeneous magnetization gives rise to the long- 
range effect provided by ±1 triplets. The relative angles 
between the spin polarizations in the superconducting 
leads and in the ferromagnet can be varied by control- 
ling the injected spin polarizations, making it possible to 
vary the magnitude and sign of the Josephson current. 
This enables manipulations of 7r-junctions. In addition 
to the Josephson supercurrent, which is driven by the 
difference in the condensate phases, there is also a dissi- 
pative DC current. The latter is induced by the spin po- 
larization flow through the ferromagnetic layer with spin 
dependent conductivity. This dissipativc current also can 
be manipulated by varying the injected polarization an- 
gles. As it will be shown, at some angles it vanishes, so 
that the dissipative and supercurrents can be measured 
independently. 

Various effects of an injected spin polarization and spin 
current on the electric transport in SFS junctions^ and 
other superconducting systems^— have been recently 
considered. Despite this interest, the fact that the ex- 
change interaction transfers the spin polarization from 
the quasiparticles to the condensate has not been ad- 
dressed so far. 

The article is organized by the following way. In Sec. II 
an expression is derived connecting the triplet compo- 
nents of the anomalous Green function to the noncqui- 
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FIG. 1: (Color online) A sketch of the system. The electric 
current flows in normal leads N through contacts with ferro- 
magnetic leads F_l and F_r. Spin density is injected from Fl 
and F_r into N and further penetrates across tunneling barri- 
ers into superconductors Si and Sr. The Josephson current 
flows between these leads through a ferromagnetic layer F. 
Arrows show possible magnetizations of the ferromagnets. 



librium spin polarization in superconducting leads. In 
Sec. Ill the Josephson and dissipative currents are calcu- 
lated. Finally, our results arc discussed in Sec. IV. 



II. TRIPLET ELECTRON PAIRING FUNCTION 
INDUCED BY SPIN INJECTION 

How to efficiently inject a spin polarization into para- 
magnetic metals is well knowni^. A noncquilibrium spin 
accumulation is induced by the electric current through 
a paramagnetic-ferromagnetic interface. We consider the 
scenario that the spin polarization further diffuses from 
a paramagnet through a resistive barrier into a supercon- 
ducting lead, so that the electric circuit where the spin 
injection takes place is effectively separated from the su- 
perconducting circuit. We assume that the steady state 
spin polarizations are generated in both superconduct- 
ing leads, in the vicinity of the F-layer. The sketch of 
the system is shown in Fig.l. For clarity, we simplify 
the problem by assuming that the FS contacts contain 
a barrier, so that the proximity effect is weak. We also 
assume that the spin relaxation time r sp i n in the leads is 
long, so that the spin diffusion length / sp ; n is large com- 
pared to the SN contact sizes and the coherence length. 
Consequently, the spin densities s^rjA and the order pa- 
rameters Al(r) only vary slowly in space near the left 
(L) and right (R) contacts. 

The electronic transport through an SFS system, 
whose characteristic dimensions are larger than the elas- 
tic mean free path, can be described in terms of Us- 
adcl equations for angular averaged Green's functions 
g (for a review seei^). These functions are matrices in 
the Keldysh, spin, and Nambu spaces. We choose the 
spin and Nambu spaces so that the one-particle destruc- 
tion operators are c lkt = c kt , C114 = C H> c 2kt 
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Nambu spinor components, while f arL d \. are the spin 
indices. The Keldysh component g K of the Green func- 
tion can be represented asi^ 



g K = g' h - hg a . 



(1) 



where g r and g a are the retarded and advanced functions, 
respectively, and the distribution function h is a diagonal 
matrix in the Nambu space. 

In order to determine the distribution h in the super- 
conducting leads, the interfaces between these leads and 
the spin-polarized normal metals must be considered. We 
use standard boundary conditions relating fluxes through 
S-N (S-F) interfaces to Green functions in superconduc- 
tors and normal metals (ferromagnets). It is assumed 
that the spin relaxation rates in the superconducting 
leads are slow enough (l sp i n r sn a s ) and the leak- 
age of the spin polarization through the SF boundary is 
sufficiently slow r sn /A sn -C r s f/A s f, where l/r sn and 
I /r s f are the interface conductances (per unit square) of 
SN and SF interfaces, A sn and A s f are the SN and SF 
contact areas, and er s is the normal-state conductivity of 
the superconductor's lead. With these assumptions, the 
distribution functions in the superconductor, h^ s \ and 
normal metal, h^ n \ are equal to each other, = hy 1 '. 
We further assume that nonequilibrium spins in N-leads 
are thermalized with chemical potentials fXf and /ij, for 
the two spin directions. Therefore, denoting by the sub- 
scripts 11 and 22 the corresponding matrix elements in 

the Nambu space, we get for 



and h 
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At the same time, the retarded (g r ) and advanced (g a ) 
Green functions have the same forms as in an equilibrium 
superconductor. 

Our calculation so far re-iterates the conventional wis- 
dom of spin-injection in superconductors: the effects are 
limited to a spin-dependent statistical distribution func- 
tion, while the retarded and advanced Green functions do 
not change. In this picture, spin injection does not lead to 
the appearance of triplet correlations in the condensate 
wave-function, which would cause long-range Josephson 
tunneling through a ferromagnetic layer. Fortunately, 
there is a mechanism to generate triplet correlations in 
spin-polarized superconducting leads, which others have 
so far overlooked. The electron-electron exchange in- 
teraction provides a coupling between a spin accumu- 
lation and the spectral properties of superconductors, 
in that spin polarized quasiparticles produce an effec- 
tive Zeeman field. The latter, in its turn, gives rise to 
triplet correlations that are described via the correspond- 
ing spin components of the anomalous functions g{ 2 and 
g1 2 . In Fermi-liquid theory, the effective Zeeman energy 
is ej; C (crN), where N is a unit vector parallel to the in- 
jected spin polarization S = NS and 



"C'-kt' where the labels 1 and 2 denote the 



GS/2N F . 



(3) 
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The spin-accumulation magnitude is 



S=- 



4(1 + G) 



diuTr 



(1 



(4) 



where r 3 and a z are the Pauli matrices acting in the 
Nambu and spin spaces, respectively, and Np is the den- 
sity of states at the Fermi level. The renormalization 
factor 1/(1 + G), where G is the exchange Landau- Fermi 
liquid parameter, appears when the spin-density of Eq. 
2] is expressed in terms of a semiclassical Green func- 
tion integrated over energy^. This factor is not quali- 
tatively important in our case, since G is not too close 
to the paramagnet instability G = — 1^ The exchange 
Coulomb interaction in metals gives rise to a negative G. 
For example, the calculated value is -0.17 in Ali&. The 
spin density ((4]) strongly depends on temperature, mostly 
via the temperature dependence of the superconducting 
gap in the energy spectrum. In order to determine S and 
A in both leads, Eq. (jU have to be solved together with 
the S-depended selfconsistency equation for A.— 

Via the effective Zeeman energy of Eq|3l the retarded 
and advanced Green functions become spin-dependent»i 
Indeed, choosing the quantization axis along N, the 

anomalous functions = g[ 2 tt an< ^ fit = — #124.4. 
come 
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where the phase <j) of the order parameter A equals (f>L 
and (j)R at the left and right contacts, respectively. The 
triplet component of this function with O-spin-projcction 
onto the z-axis is /q = (fL + /L.)/\/2, while the triplet 
components with ±l-projections vanish, f±i = f^nn = 
0. The advanced function, as well as the conjugated func- 
tions f\ are determined from symmetry relations. 

It is more transparent to discuss the Green functions in 
a basis where the spin quantization axis is parallel to the 
magnetization in the ferromagnetic layer, which is along 
z, as shown in Fig. 1. We assume that the spin polariza- 
tions in the left and right leads are rotated with respect to 
this axis by the angles 9l and 9 R , respectively. We follow 
the convention that the three components of the triplet 
/0)/ii f-i are related to a 3D vector a = (a Xl a y ,a z ) with 
a z = fo, a x = (/_! - h)/V2 and a y = + fi)/V&L. 
Hence, in the geometry shown in Fig. 1, after a rota- 
tion of a around the j/-axis, we get in the new basis 
fi = Jocose and f[ = = -f sm9/V2. So, by 

using Eq. ([5]) the triplet components in the left and right 
superconducting leads arc 



f±lR(L) 



sin 9 



R(L) 



(In + fit) : 



(0) 



where the labels r and a have been omitted from here and 
the same magnitudes of e xc are assumed in both leads. 



In the new basis, the distribution function (|2|) is 

, , {l + a z cos8 L(R) ) {1- <t z cos6 l{r) ) 
n-LiR) = «t ~ r ti_\ - 



a y sm9 LiR) ) 
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III. 



THE JOSEPHSON AND DISSIPATIVE 
CURRENTS 



What we have established is that the superconduct- 
ing leads acquire triplet pairing correlations determined 
by non-equilibrium spin polarizations whose directions 
are tilted with respect to the ferromagnet's magnetiza- 
tion in the SFS junction. We will show that the current 
through such a triplet pairing-ferromagnet-triplet pairing 
system consists of two parts, a dissipative contribution 
controlled by the non-equilibrium distribution of spins 
in the device, and a super-current driven by the phase 
difference between superconductors and provided by the 
triplet components of the superconducting condensates 
in the left and right leads. 

Let us first consider the dissipative current. It can be 
expressed in terms of the distribution function hf inside 
the ferromagnet. Due to precession in the exchange field 
B ox , the spins that are not parallel to it decay quickly 
on the length-scale ^/Z?f/B cx , where Df is the diffusion 
constant. Therefore, only the components of h{ that are 
parallel and anti-parallel to z, denoted as hff and hfi, re- 
main finite inside the ferromagnet, if the junction length 
L ^> sj D[/B CK . When the spin relaxation length is 
larger than L, in the linear approximation these collinear 
components satisfy the spin-conserving diffusion equa- 
tion Df a W x hf a = 0, where a =t,|, that takes into ac- 
count spin-dependent diffusion coefficient in a strong fer- 
romagnet. The solution of this equation is a linear func- 
tion of x whose slope is obtained from the boundary con- 
ditions T^sfcrCTfcrV ' x hf<j\ x = XL{R) = /lL(_R)o- - h ta \ x=XL{R) , 

where h^^fj^ are given by the first two terms of Eq. (JT]). 
Taking into account that r s { a and af a can depend on the 
electron spin and assuming equal barrier transmittances 
at L and R contacts we obtain 



1 + 2 7ct L' 



(8) 



where Xr(l) = ±L/2 and j a = (r s f a a f a / L) > 1. Using 
Eqs. (HJ and (JT]) we compute the dissipative part of the 
current through the junction: 



Jd 



5fi 



J duj<T f(J \7 x hf a = 



(Tff 



0-f4. 



eL\l + 2^ l + 2 7i 



(cos 9 R — cos 9l) ■ (9) 



This current is proportional to the difference in the spin- 
up and spin-down conductances (2r s f a + L/crfg.) -1 of the 
total ferromagnetic layer, including the interfaces; this 
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is the well knowni 9 - connection between spin and elec- 
tric transport in ferromagnets. The electric current at- 
tains its maximum when cos 9r — — cos 8l = il, and 
vanishes at Or = 6l, as well as at 0r,0l = ±7r/2. 
Such an angular dependence has a simple physical ex- 
planation. The electric current ([9]) is proportional to 
the spin-current through the junction. The latter attains 
its maximum when the nonequilibrium spin polarizations 
in the ferromagnetic layers are oppositely directed and 
it vanishes if these polarizations are collinear and have 
equal magnitudes. The spin current obviously also van- 
ishes if these polarizations are perpendicular to the ferro- 
magnetic magnetization axis, since perpendicular compo- 
nents do not penetrate deep into ferromagnet. The spin 
flow through the junction is accompanied by energy dis- 
sipation. It is determined by the Ohmic losses in the 
ferromagnet during transport of spin polarized electrons 
between the leads having spin dependent electrochemical 
potentials. The dissipative current of Eq. ^ is inde- 
pendent of the superconducting phases 4>l and <Pr. We 
assume that the electric potentials of both contacts are 
equal. If the load is present in the circuit, the spin current 
will induce a voltage difference. The latter, in its turn, 
can cause periodic oscillations of the Joscphson current^. 

When yj Df/B ex is much shorter than the junction 
length L and the coherence length, the up and down- 
spin Fermi surfaces become decoupled. In this regime, 
the supercurrent j s through the junction is determined 
by the decoupled tunneling of ±1 triplet Cooper pairs at 
their respective ferromagnet 's Fermi surfaces. Unlike the 
dissipative current, the spin-dependence of the electron 
diffusion coefficients and conductivities is not so impor- 
tant, at least in the case when the i? GX <C Ep. Therefore, 
in the leading approximation we set Dff = = Df, 
and a similar relation for the conductivities. Further- 
more, in the linear approximation, only the first term of 
Eq. ((HJ has to be taken into account. Moreover, since 
the Josephson current is determined by part of the dis- 
tribution function that is odd in frequency, from Eqs. 
©, © and © only the spin-independent part htf + /iq, 
contributes to the current. It is given by 
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(/iff + hfi) 



(10) 



where (m = ±1) are the retarded triplet compo- 
nents of the anomalous function in ferromagnet and 
fiM = while /£( W ) = /£(- w ). When the 

spin- relaxation length is much larger than L, and within 
the linearized approximation, obey^ 



DfVlf m + 2kjf m = , 



(11) 



with the boundary conditions^ r s f<jf\7 x f m \ x=XR{L) = 
±f m R(L), where f mR ( L ) are given by Eqs. flS©. Af- 
ter transforming the integral in Eq. (|10[) into a sum over 



I/In 



r=0.39 T 



r=0.56T 



T=0.73 T 





FIG. 2: (Color online) The Josephson current as a function 
of the spin potential S/j, measured in units of the unperturbed 
superconducting gap Ao, at L = £, where £ = vADf/2|Ao| 
and I = (2(nk B T/er%a f )xW- 2 



the frequencies w„ = TrksTftn + 1), j s can be finally 
represented in the form 

j s = sm(0 £ - 4> R ) —- 2 K , (12) 



where 



X=|A|2fcsT S kv L^h(KL) 

UJn>0,V — ±1 



(13) 



and k v = y // 2(oj n + iis5(i)/Df, with H u = e xc + v8(i and 
26 (i = (if - Hi- 



IV. DISCUSSION 



As follows from Eq. (|T2j) . the Josephson current de- 
pends on the directions of the nonequilibrium spin po- 
larizations in the superconducting leads. The current 
reaches its maximum when the spin accumulations in the 
leads are perpendicular to the magnetization in the fer- 
romagnet, Or = 9l = 7r/2. It reverses its sign when 
the spin polarization in one of the leads flips its direc- 
tion. Therefore, in the setup shown in Fig. 1, the junc- 
tion can be switched into the 7r-state by simply revers- 
ing the electric current through one of the FN contacts. 
It should be noted that the dissipative current given by 
Eq.© vanishes when the relative angles are such that 
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the supercurrcnt reaches its maximum. Hence, the dissi- 
pative transport can be turned off, a feature that can be 
important for practical purposes. Eqs. (fT2]) and (fT3|) also 
imply that the long-range proximity effect, described via 
j s , vanishes when the exchange interaction e xc = 0. The 
dependence of j s on the spin-potential 8[i is shown in Fig. 
2. A finite spin-potential causes variations of the order 
parameter A and spin density entering in Eq. (|12[) which 
have been found from a pair of self-consistent equations. 
In our calculation of A and S, we neglected the exchange 
field e xc , assuming that e xc -C <5/i. This is a realistic 
assumption, taking into account that S < SfxNp in Eq. 
([3]) and |G| is considerably less than 1 in some supercon- 
ducting metals (e.g. Al). In this limit, the dependence of 
A on 5[i is formally the same as in a thermally equilib- 
rium superconductor subject to a Zecman splitting equal 
to Sfi. Such a scenario is well studied in the literature 
(see e.gj2i). In Fig. 2, we see that the critical current 
changes sign at some values of Sfi. This is caused by 
injection of nonequilibrium spins into the ferromagnetic 
layer. As a result, the distribution function in Eq. [TTJ1 is 
different from the equilibrium distribution. There is some 
similarity of this effect with a current reversal observed 
in Joscphson transistors^.. At the lower temperature, 
the supercurrent versus spin-potential is more peaked in 
the range of higher Sfi, since the spin density increases 
sharply together with e xc in this range. At even higher 
Sfi, the superconductivity is destroyed by spin injection. 
That causes a sudden drop of the current. We believe 



that this narrow range can be easily observed experi- 
mentally in the set up shown in Fig. 1, because Sfi can 
be fine tuned by varying the current through the normal 
leads. 

Fig. 2 is calculated at L = y / D{/2\A \, that is the 
characteristic length of the ±l-triplet proximity effect 
in the range of temperatures considered. This length is 
obviously much larger than the s-wave Cooper pair pen- 
etration depth y/ Df/B cx and therefore clearly demon- 
strates how the range of the proximity effect becomes 
much longer by spin injection into the superconducting 
leads. 

In conclusion, spin injection into s-wave superconduc- 
tors can dramatically increase the stationary Joscphson 
current in SFS system. This enhancement is provided 
by ±1 triplet components of the electron pairing func- 
tion. They are generated in superconducting leads by 
exchange fields that are noncollinear with the ferromag- 
net magnetization. These fields, in turn, are induced by 
an injected spin polarization. Besides a strong effect on 
the Josephson current, spin injection also gives rise to a 
dissipative current that at zero bias potential is induced 
due to spin dependence of the ferromagnct conductivity. 
Both Josephson and dissipative currents can be manipu- 
lated by varying the injected spin directions in the leads 
enabling control of 7r-j unctions. 

Anatoly Mal'shukov gratefully acknowledges the hos- 
pitality of NTNU. 
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